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Abstract 

Possible rotation states of two satellites of Saturn, Prometheus (S16) and Pan- 
dora (S17), are studied by means of numerical experiments. The attitude stability 
of all possible modes of synchronous rotation and the motion close to these modes is 
analyzed by means of computation of the Lyapunov spectra of the motion. The stabil- 
ity analysis confirms that the rotation of Prometheus and Pandora might be chaotic, 
though the possibility of regular behaviour is not excluded. For the both satellites, the 
attitude instability zones form series of concentric belts enclosing the main synchronous 
resonance center in the phase space sections. A hypothesis is put forward that these 
belts might form "barriers" for capturing the satellites in synchronous rotation. The 
satellites in chaotic rotation can mimic ordinary regular synchronous behaviour: they 
preserve preferred orientation for long periods of time, the largest axis of satellite's 
figure being directed approximately towards Saturn. 

1 Introduction 



In 1980s, IWisdom et al.l (119841 ) and IWisdoml ( 119871 ) demonstrated theoretically that a plan- 
etary satellite of non-spherical shape in an elliptic orbit can rotate in a chaotic, unpre- 
dictable way. They found that the most probable candidate for the chaotic rotation, due 
to pronounced shape asymmetry and significant orbital eccentricity, was the satellite of 
Sat urn Hyper i on (S7 ) Later on, a direct modelling of observed light curves of H yper- 
ion dKlavetterl . ll989al Jbl: iThomas et all Il99,4 black et all Il99.4 iDevyatkin et all l2002h con- 
firmed the chaotic character of its rotation. 



Recently it was found in a theoretical research (jKouprianov and Shevchenko! . 120051 ) . that 
two other satellites of Saturn, Prometheus (S16) and Pandora (S17), can also reside in a 
state of chaotic rotation. Contrary to the case of Hyperion, chaos in rotation of these two 
satellites is due to fine-tuning of the dynamical and physical parameters rather than simply 
to a large extent of a chaotic zone in the rotational phase space. 

It is remarkable that the orbital dynamics of Prometheus and Pandora are known to be 
chaotic wit h the Lyapunov time of only three y ears, and this dynamical chaos is directly ob- 
serve d; see ( IGoldreich and Rappaportl . l2003al ibl; ICooper and Murray! . 120041 : Farmer and Goldreich 
20061 ) . However note that the theoretical inferences on chaos in rotation of these satellites 
are completely independent from the existence of orbital chaos. 
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In the present paper we study the problem of rotational dynamics of these two satel- 
lites in detail, exploring the attitude stability not only in the centers of synchronous reso- 
nances, but also in the phase space in the vicinities of the resonances. Our analysis includes 
period-doubling bifurcation modes of synchronous spin-orbit states. We take into account 
all available modern observational data. 

An important problem, first of all from the observational point of view, is the follow- 
ing: whether there exists a preferred orientation of the satellites rotating chaotically, or all 
orientations are equiprobable? The orientation of a satellite in chaotic rotation, generally 
sp eaking, is not n ecessarily isotropic. This was demonstrated in a numerical experiment 
by IWisdoml (119871 ) in calculating the rotatio n of Phobos in t he vicinities of the separatrices 
of the 1:2 spin-orbit resonance; see Fig. 5 in ( IWisdoml . 119871 ). 

The plan of the present work is as follows. First, we study the attitude stability of 
the planar rotation of Prometheus and Pandora with respect to tilting the axis of rotation. 
The stability analysis is carried out for all possible exact modes of synchronous resonance 
(a-resonance, (3- resonance, and the period- doubling bifurcation mode of a-resonance), as 
well as for the trajectories of all possible kinds (periodic, quasiperiodic, chaotic in the planar 
problem) on a representative section of the phase space of planar rotation. Then the problem 
on the preferred orientation of Prometheus and Pandora in chaotic rotation is considered. 



2 The reference frame and the equations of motion 



We suppose that a satellite represents a non-spherical rigid body moving around a planet 
(a gravitating point) in a non-perturbed elliptic orbit with the eccentricity e. Location of 
the satellite in the orbit is determined by the true anomaly / or the eccentric anomaly E. 
The shape of the satellite is described by a triaxial ellipsoid with the principal semiaxes 
a > b > c and the corresponding principal central moments of inertia A < B < C . The 
size of the satellite is much less than the radius of the orbital motion. The dynamics of the 
three-dimensional rotation of the satellite is determined by the parameters e, b/a, c/b and 
the initial conditions of the motion. The angular velocities of rotation are expressed in the 
units of the orbital mean motion, and the "satellite — planet" distance is expressed in the 
units of the semimajor axis of the orbit. One orbital period corresponds to 27r time units. 

The oblateness of Saturn causes precession of the pericenters of the orbits of the satellites, 
the precession rate being equal to 3. 1911 x 10~ 5 °s~ 1 in the case of Pro metheus and 3.0082 x 



10 -5 in the case of Pandora (jGoldreich and Rappaportl . l2003bl ). Hence the periods 
of precession are equal to 0.36 and 0.38 yr, respectively. This is much greater that the 
Lyapunov times (less than 1 d; see Kouprianov and Shevchenko f 2005h ) of the rotation 
of these satellites, if it were chaotic. Since the timescales are so different, we ignore the 
precession of orbits in our study. 

For describing the orienta t ion of the satellite we use an inertial f rame identical to that 
used in (IWisdom et all I1984J : iMelnikov and Shevchenko! . Il998l . 120001 ) . This Oxyz frame is 
defined initially at the pericenter of the orbit as follows: the x axis is directed along the 
"orbit pericenter — planet" vector, the y axis is parallel to the vector of the orbital velocity 
at the pericenter, the z axis is orthogonal to the orbital plane and completes the reference 
system to a right-handed system. Orientation of the satellite with respect to the axes of 
the Oxyz frame is defined by a sequence of imaginary rotations of the satellite by the Euler 
angles 9, <f), ip from an initial position until the satellite reaches its actual orientation. In 
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the initial position the axes a, b, c coincide with the axes x, y, z, respectively. The axes a, 
b, c are directed along the principal axes of inertia with the moments A, B, C, respectively, 
and are "frozen" in the satellite. The imaginary rotations are carried out in the following 
sequence: first, rotation by 9 about c, second, rotation by about a, and third, rotation by 
■0 about b. 

The definit i on of the Euler angles adopted here is identical to that described and used by 
Wisdom et al.l (119841 ) and is different from the usual one. The reason for using the alternative 
Euler angle set is that the standard one has a coordinate singularity at the point where a 
satellite's c axis (the axis of the maximum moment of inertia) is orthogonal to the orbit 
plane. In the adopted frame this position corresponds to <fi = 0, while the singularity is 
shifted to <fi = ±7r/2. The latter value corresponds to the satellite's axis of rotation lying in 
the orbit plane. 

Rotation of the satellite i s described by Euler 's dyna mic and kinematic equations. The 
dynamic equations ( Beletsky . 1965 ; Wisdom et al. . Il984"l ) can be written as 
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Here G is the gravitational constant; M is the mass of the planet; u a , u;&, u c are the pro- 
jections of the vector of the angular velocity cJ on the axes a, b, c; r = a(l — ecosE) is the 
"satellite — planet" distance, a is the semimajor axis of the orbit; a, (3, 7 are the direction 
cosines of the principal axes of inertia with respect to the direction to the planet. 

The kinematic equa tions and expressions for direction cosines, in the re ference frame 
used here, according to (IWisdom et al.l . Il984 ; iMelnikov and Shevchenkol . ll998l ). are given by 
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Eqs. ([I]), (T5]), and ([3]) are used in what follows for calculation of the rotational dynamics. 



3 



3 Synchronous resonance regimes and attitude stabil- 
ity of rotation 

3.1 The "Amalthea effect" 



In (IMelnikov and Shevchenkol . Il998l . 120001 ) it was found for Amalthea (J5) that two differ- 
ent synchronous regimes of rotation, "q-resonance" and "/3-resonance" , coexist in the phase 
space of rotation of t his satellite. By means of calcu lation of the Lyapunov characteris- 



tic exponents (LCEs), iMefnikov and Shevchenkol (119981 ) showed that the planar rotation of 



Amalthea in the /3-resonance center and its neighbourhood was stable with respect to tilting 
the axis of rotation, while its rotation in the q-resonance center and its neighbourhood was 
unstable. (Note t hat the terms "q-resonanc e " and "/9-resonance" had not yet been used 



in that paper.) In (IMelnikov and Shevchenkol . |2000| ) by means of calculation and statistical 
analysis of the multipliers of the periodic solutions corresponding to the a- and /3-resonances, 
the stability of rotation of Amalthea at the centers of the both modes was investigated on 
a grid of values of the inertial parameters. The conclusion was made that Amalthea could 
not reside in q-resonance. 

Let us introduce the parameter ojq = 



3(B — A)/C, which is the frequency of small 
amplitude oscillat ions of a satellite in synchronous resonance (see ( iWisdom et all [1984; 



Shevchenkol . Il99flh ). It roughly characterizes the dynamical asymmetry of the satellite shape. 
If we accept the data of lEA'99l . then w = 1-058 for Prometheus and uq = 0.812 for Pandora. 

For two centers of synchronous resonance to coexist in the phase space of rotational 
motion (t he "Amalthea effect"), the ojn p arameter must exceed unity slightly; see Figs. 1 
and 2 in ( IMelnikov and Shevchenkol . |2000| ) and also Fig. El in the present paper. Note that, 
in addition, the orbital eccentricity must obey a certain constraint: it must not be too high. 
The limit on the eccentricity, which allows coexistence of the periodic solutions corres ponding 
to q-resonance and /3-resonance, is as follows: e < ^-(u — l) 3 / 2 ( Beletsky ( 1965 . Ch. 2); 
Markeevl (Il990l p. 366)). 

The "Amalthea effect" takes place for prolate satellites. Indeed, cuq = y3(B — A)/C = 



3 (a 2 — b 2 )/ (a 2 + b 2 ) for a triaxial ellipsoid with homogeneous density; see ( iKouprianov and Shevchenk 

20061 . p. 396). The parameter ujq > 1 if c < b < a/y/2. In other words, two semiaxes should 
be less than 0.7 of the third one. The "Ama l thea effect" was considere d and discussed 
in detail in ( Kouprianov and Shevchenko . 20061 : IMelnikov and Shevchenkol . 2007 ). It was 
shown that this effect might be abundant amongst minor planetary satellites (the satellites 
with diameters less than 100 km) moving in close-to-circular orbits. 



3.2 Synchronous states: a-resonance, /3-resonance and period- 
doubling bifurcation mode of a-resonance 

The available data on the parameters of the figures of Prometheus and Pandora are collected 
i n Tab le HI Note that the data on the sh apes of Prometheu s and Pandora derived bvlThomas 
( 1989) are tabulated in the reports by Seidelmann et al. ( 2002 ). Seidelmann et al.T "( 2005 ). 



and ISeidelmann et al.l (120071 ). Besides, the data on the adopted positions of the poles of 
Prometheus and Pandora can be found in these reports; these data are the same in the three 
references. We set the or bital eccentric ity equal to e = 0.002 for Prometheus and e = 0.004 



for Pandora, as in lEA'99l . According to lGoldreich and Rappaportl ( I2003al Figs. 5 and 6), on 
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Table 1: The shape parameters for Prometheus and Pandora 



Prometheus (S16) Pandora (S17) References 



b/a 


c/b 


b/a 


c/b 




0.714 


0.740 


0.764 


0.786 


Wisdom r 19871 


0.676 


0.680 


0.800 


0.705 


Thomas ( 1989) 


0.586 


0.706 


0.737 


0.738 


Stooke ( 19931 


0.608 


0.726 


0.741 


0.782 


Gozdziewski and Macieiewski (1995) 


0.676 


0.680 


0.800 


0.705 


EA'99 


0.734 


0.696 


0.773 


0.804 


Porco et al. (2006) 



the timescale of 20 years the eccentricity of Prometheus varies in the limits 2.27 x 10 3 
2.30 x 10~ 3 , and that of Pandora in the limits 4.35 x 10~ 3 - 4.38 x 10" 3 ; see also (French et al. 



2003 ). 

The "e — co " diagra m for the satellites with known valu es of the shape parameters, based 
on the data compiled in ( IKouprianov and Shevchenkd . l2005l ). is shown in Fig. Ufa). The theo- 
retical boundaries of the zones of existence of a-resonance, /3-resonance, an d period-doublin g 
bifurcation mode aut of a- resonance are indicated according to the data in (jMelnikovl . l200ll ) . 

As follows from Fig. QJa), period-doubling bifurcation mode a^a of a- resonance can 
be present in the phase space of rotation of Pr ometheus and Pandora. For Pandora, this 
confirms an earlier analysis by iMelnikovl ( 120011 ) made on the basis of a single estimate of 
the inertial parameters. In Fig. QJb), a part of the "e — ojq" diagram with location of 
Prometheus and Pandora indicated according to the data of Table [1] is given in higher 
resolution. In the case of Pandora, /3-resonance does not exist for all the data co nsidered. 
In the case of Prometheus, /3 -resonance does not exist only for the data due to IWisdom 
(119871 ) and iPorco et al.l (120061 ). According to all other sources, Prometheus lies in the zone 
of existence of this mode. For all the data on Prometheus and Pandora, period-doubling 
bifurcation mode of a-resonance exists. 

In Fig. |2l the phase space section of the planar rotational motion of Prometheus is shown. 
The section is defined at the pericenter of the orbit; i.e., the variables are mapped each orbital 
period. Note that the planar problem (that with <f) = = 0) has one and a half degrees of 
freedom. The center of a-resonance (the lower one in the section) and that of /3-resonance 
(the upper one) are indicated in the section. Besides, there is period-doubling bifurcation 
mode «bif, that manifests itself in the two prominent regular islands inside the chaotic layer 
at the left and at the right sides from the center of /3-resonance. 

We see that Pr ometheus i s subj ect to the "Amalthea effect" . In a different terminology 



this was noted in (jMelnikovl . 1200 ll ). What is more, period-doubling bifurcation mode of 
a- resonance is possible for this satellite, i.e., the potentially possible rotational dynamics are 
very rich. 



3.3 Attitude stability of exact synchronous rotation 

Let us consider the stability of synchronous rotation of Prometheus and Pandora with respect 

to tilting the axis of rotation i n the cases of exact a-resqnanc e, /^- resonance, and period- 

doubl ing bifurcation mode a hii . iMelnikov and Shevchenkol (l2000h and IKouprianov and Shevchenkol 
(120051 ) performed a research of the stability on the {A/C, B/C) plane. Here we use another 
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Figure 1: (a) Location of the satellites with known values of the shape parameters in the 
"e — u " diagram, (b) Location of Prometheus (crosses) and Pandora (circles) in the 
diagram "e — w " according to Table [TJ 
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Figure 2: The phase space section of the planar rotational motion of Prometheus. The exact 
synchronous resonance modes are indicated. 

plane, namely, the (c/b, b/a) planeQ, which is more graphical. 

In the cases of a-resonance and /3-resonance we use a method based on the analysis 
of the modal structure of the differential distribution of the computed modules of multi- 
pli ers of periodic solutions of the e quations of motion. The method is described in detail 
in ( Melnikov and Shevchenko . 2000l ). 

In the considered problem, the system of equations of motion in variations with respect to 
the periodic solution consists of six linear differential equations of the first order with periodic 
coefficients. Numerical integration of the syst em allows one to obt ain the matrix of linear 
transformation of variations for one period; see (IWisdom et al.l . ll984l ). The periodic solutions 
in the given problem are characterized by three pairs of multipliers. The distributions of 
the modules of multipliers are built for a set of trajectories corresponding to a center of 
synchronous reso nance on a grid of values o f the b/a and c/b parameters. Analysis of the 



distributions (see iMelnikov and Shevchenko! (120001 )) allows one to separate orbits stable with 



respect to tilting the axis of rotation from those which are unstable. 

Analysis of the attitude stability of period-doubling bifurcation mode a^i is carried out 
by means of computation of the whole spectrum of the Lyapunov characteristic exponents 
(LCEs) for a set of values of the b/a and c/b parameters, followed by analysis of the differ- 
ential distributions of the computed values of the LCEs. By means of a similar method we 
investigated the stability of rotation of planetary satellites (IMelnikov and Shevchenko! . Il998l ) 
on sets of initial data of the trajectories. Here the distributions of the LCE values for the 
trajectories with the initial data taken at exact period-doubling bifurcation mode «bif are 
built on a grid of values of the b/a and c/b parameters. We set the grid resolution equal to 
0.001 in both b/a and c/b axes. Analysis of the modal structure of the distributions allows 
one to separate the stable trajectories, for which all three indices are zero, and the unstable 



: We are grateful to A. Dobrovolskis for the advice to use the given coordinates. 
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Figure 3: The d8/dt coordinates of the centers of synchronous resonance in the phase space 
section in dependence on the u parameter. The curves starting on the left side of the plot 
correspond to a-resonance, those starting on the right side correspond to /3-resonance. 



ones for which at least one of the LCEs is distinct from zero. 

The true values of the LCEs are supposed to be the limits of the computed values when 
the time of computation tends to infinity. The time of computation is necessarily finite. 
However, it is implied henceforth that the obtained numerical values in the case of chaotic 
trajectories represent the true LCE values, because the time of computation was taken to be 
long enough for the computed LCEs of the chaotic trajectories to saturate (i.e., increasing 
the computation time would not make the computed LCE values less; a "plateau" is reached 
in each case). In what concerns the regular trajectories, the obtained numerical values of 
the LCEs tend to zero with increasing the computation time. 



For the computation of the LCE spectra , we use the HQRB method (yon Bremen et al. 



1997), programmed as a software package in (IShevchenko and Kouprianov 



2002 



Kouprianov and Shevc 



2003). The Dormand-Prince integrator DOPRI8 ( Hairer et al. . 19871 ). realizing in Fortran 



the 8th order Runge-Kutta method with the step size control, is used for integration of 
equations of motion (JTJ) , (J2j) . 

In Fig. |3j the d8/dt coordinates of the centers of synchronous resonance in the phase 
space section in dependence on the uq parameter are shown for the values of the orbital 
eccentricities of both satellites. The curves starting on the left side of the plot correspond to 
a-resonance, those starting on the right side correspond to /3-resonance. As it is clear from 
Fig. EJ a- and /3-resonances coexist in a substantial interval of cuo, if the orbital eccentricities 
are so small. /3-resonance is born at Uq ~ 1; a- resonance disappears at a large value of ujq, out 
of the presented plot limit. The problem of bifur cations causing the birth a n d disa ppearance 
of a- and /3-resonances is considered in brief by iMelnikov and Shevchenkol fl2000l ). 

The computed regions of stability and instability are shown in Fig. [J] (for Prometheus) 
and in Fig. (for Pandora). The regions of stability are shown in light gray, the regions 
of minimum (one degree of freedom) instability are shown in dark gray, and the regions of 
maximum (two degrees of freedom) instability are shown in black. The lines of constant 
value of the Uq parameter are depicted for orientation. 

It follows from the diagrams in Figs. Ufa) and (b) that a-resonance is attitude unsta- 
ble for most of the observational data for Prometheus . Earlier this instability was noted 
in ( Melnikov . 2001 : Kouprianov and Shevchenkol . 2005 ) on the basis of a single estimate of 
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Figure 4: Regions of stability and instability with respect to tilting the axis of rotation; e 
0.002 (Prometheus): (a) for the center of a-resonance, (b) for the center of /3-resonance, (c) 
for the exact period-doubling bifurcation mode of a-resonance. The locations of Prometheus 
according to the data of Table [1] are indicated by crosses. 
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0.0-| , 1 , 1 , 1 , 1 , 

0.0 0.2 0.4 0.6 0.8 1.0 

clb 

b) 

Figure 5: Regions of stability and instability with respect to tilting the axis of rotation; 
e = 0.004 (Pandora): (a) for the center of a-resonance, (b) for the exact period-doubling 
bifurcation mode of a-resonance. The locations of Pandora according to the data of Table Q] 
are indicated by crosses. 
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the inertial parameters. In the case of the data due to iPorco et al.l ( 120061 ). a- resonance is 
close to instability. Rotatio n of Prometheus in /3-resonan c e is c lose to instability, in agree- 
ment with an inference by iKouprianov and Shevchenkol (120051 ). In the case of Pandora 



(see Fig. Ef a)) , rota ti on in q-resonance is close to instability for the dat a of (| Wisd om 



1987; Thomas, Il989l: IPorco et all l200fih and is unstable for the data of (IStookeL 119 93 



Gozdziewski and Maciejewskil . ll995l ). in general agreement with (jKouprianov and Shevchenkol . 
2005h . where the stability analysis was based on a single estimate of the inertial parameters. 

From the diagrams in Fig. 0](c) and Fig. E3(b) it is clear that synchronous rotation in the 
«bif mode is attitude unstable for both Prometheus and Pandora. 



3.4 Attitude stability in the general case 

Not only the stability in exact synchronous rotation is of interest. Let us study the attitude 
stability of the trajectories in the vicinities of the exact synchronous states on representative 
sets of initial data. Computation of the LCEs on a grid of initial data, followed by analysis of 
the distributions of the computed LCEc, enables one to accompl ish such a study. A similar 
study with the use of the maximum LCEs was carried out in ( jMelnikov and Shevchenkol . 
19981 ) for Phobos, Deimos, Amalthea, and Hyperion. Here we study Prometheus and Pan- 
dora by means of analysis of the LCE spectra. Besides, we significantly increase the resolution 
of the initial data grid, as well as the time interval on which the LCEs are computed. For 
the surface of section we choose the (9, d9/dt) plane taken at t — 2Trm, m = 0, 1, 2, . . ., i.e., 
defined at the o rbit pericenter. In the computations we adopt the values of b/a and c/b as 
given in 



EA'99 



(see Table 1). 

The computation of the LCEs has been carried out for two sets of trajectories, (i) and 
(ii), defined by the following choice of initial data: (i) 9 = 0, (ii) 9 = 7r/2; and in the both 
sets d9/dt is taken in the range from —0.5 to 2.5 with the step equal to 0.003 in the case 
of Prometheus, and in the range from 0.0 to 2.0 with the step equal to 0.002 in the case of 
Pandora. The initial conditions also include <p = ip = 0, d<p/dt = dip/dt = 0; the motion 
starts (i.e., t = 0) at the pericenter. Thus the LCEs have been computed for 2000 trajectories 
for each satellite. 



Following iMelnikov and Shevchenkol (119981 ) . in order to separate the regular and chaotic 
orbits, we build differential distributions of the computed LCE values. The distribution has 
two peaks; one of them corresponds to the chaotic orbits, and one to the regular orbits. 
On increasing the time interval of integration, the peak corresponding to the chaotic orbits 
remains motionless, while the abscissa of the peak corresponding to the regular orbits tends 
to zero (or to minus infinity in the logarithmic scale). Thus the sets of regular and chaotic 
trajectories are separated. The abscissa of a point between the peaks gives the numeric 
criterion for separation of the sets. Increasing the time interval of integration allows one to 
make this numeric criterion more precise. 

The LCE dependences on the initial value of d9/dt in the case of 9 = are shown 
in Fig. [6] They have been computed on the time intervals t = 10 4 and 10 5 . From the 
plots it is clear that with increasing the integration time the computed LCE values for 
the chaotic trajectories remain constant, while the computed LCE values for the regular 
trajectories decrease. The horizontal dashed lines in the plots correspond to the adopted 
value of the LCE decimal logarithm (equal to —3 in all cases) separating the chaotic and 
regular trajectories. This criterion has been derived by means of building the distributions 
of the computed values of the LCEs at t = 10 5 . It separates the peaks corresponding to the 
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Figure 6: The dependence of the LCEs on initial data, (a) e = 0.002, uq = 1.058 
(Prometheus); (b) e = 0.004, uq = 0.812 (Pandora). The gray curves correspond to the 
computation time t = 10 4 , the black curves to t — 10 5 . The horizontal dashed lines corre- 
spond to the adopted value of the LCE decimal logarithm (equal to —3 in all cases) separating 
the chaotic and regular trajectories. 
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chaotic and regular trajectories in the distributions. 

The separation of the regular and chaotic orbits, accomplished by means of analysis of 
the LCE distributions, allows one to build two variants of the phase space section, one for 
all trajectories, and the other one only for the attitude stable trajectories. These double 
variants are shown for Prometheus in Figs. E(a) and^b), and for Pandora in Figs. [S^a) 
and Mb). The phase space sections for all trajectories are shown in Figs. [TJ^a) and IHa), 
while those for the attitude stable trajectories solely are shown in Figs. [7(b) and[8]^b). In 
order that structure of the section containing all trajectories were clearly discernible, it is 
constructed with a relatively low resolution of the initial data grid. 

In the case of Prometheus (Fig. [21 Fig. Ella)), there are two centers of synchronous 
resonance: a-resonance (the lower one in the section) and /3-resonance (the upper one). 
Besides, there exists period-doubling bifurcation mode at,if, located inside the chaotic layer 
at the left and right sides of the center of /3-resonance. The librational trajectories are 
attitude stable only in the nearest neighbourhood of the center of /3-resonance. Alternation 
of ring-like zones of stable and unstable motion is clearly seen for the librational trajectories 
enclosing the major synchronous state. This alternation is also clearly noticeable in Fig. |6^a), 
where the computed values of the LCEs are given in function of the initial d6/dt value. A 
broad band of quasiperiodic trajectories under the lower branch of the basic chaotic layer is 
attitude unstable. On the contrary, the motion above the upper branch of the layer is stable. 

In the case of Pandora (Fig. |8]), there exists period-doubling bifurcation mode abif located 
inside the chaotic layer at the left and right sides of the center of a-resonance. Alternation 
of ring-like zones of attitude stable and unstable librational trajectories enclosing the main 
synchronous state is even more pronounced than in the case of Prometheus. 

During the process of tidal capture in synchronous resonance, both satellites inevitably 
cross these intermittent belts of attitude instability. When such belts are present, attaining 
exact resonance might be more difficult in comparison with the usual situation when they are 
absent, because the satellite would tend to deviate from planar rotation in these zones, due 
to the attitude instability. So, one can put forward a hypothesis that these belts might form 
"barriers" for capturing the satellites in synchronous rotation. Some numerical-experimental 
as well as theoretical work is necessary to infer whether this hypothesis is right or not. In 
particular, it is necessary to compare the timescale of developing the attitude instability in 
a belt with that of crossing the belt due to tidal evolution. 



4 Preferred orientation in chaotic rotation 



In a theoretical research flKouprianov and Shevchenkol . l2005h it was found that Prometheus 



and Pandora are likely to be in a state of chaotic rotation. An important problem is whether 
there exists a preferred orientation of the satellites in chaotic rotation, or their "chaotic 
tumbling" is isotropic? This is important for drawing conclusions about the character of 
rotation from observational data. 



For describing rotation of a satellite we use a set of Euler angles adopted in (jWisdoml . 



19871 ) . It is different from that used above. The reason for the change is that the anisotropy 
of orientation with respect to the direction to the planet is described straightforwardly in 
the new set. The difference between the old set and the new one consists in the sequence 
of imaginary rotations by the angles from the initial position (identical to that in the old 
system; see above) to the actual orientation of the satellite. In the new set, the rotation is 
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Figure 8: The phase space section for e = 0.004, u = 0.812 (Pandora): (a) all trajectories, 
(b) only attitude stable ones. 
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made first by 6 about c, second, by <p about b, third, by — ip about a, until the axes of inertia 
of the satellite coincide with the actual orientation. 

Therefore the angle <fi in the new set is the angle between the largest axis of satellite's 
figure (the axis of the minimum moment of inertia) and the orbit plane, and the angle (9 — f) 
is the angle between the direction to the planet and the plane containing the largest axis of 
sat ellite's figure an d orthogonal to the orbit plane. 



Wisdom! (119871 ) constructed projection of a chaotic trajectory of spatial rotational motion 
of Phobo s to t he plane (0, 9 — nt), where n is the orbital mean motion (see Fig. 5 in 
(jWisdoml . 119871 ); note that 9 — f « 9 — nt for small eccentricities). The rotation of Phobos 
with the model initial conditions close to the separatrices of the 1:2 spin-orbit resonance 
was considered. Planar rotation of Phobos is unstable with respect to tilting the axis of 
rotation not only near the separatrices of this re s onanc e, but also p ractically in th e whole 
1:2 resonance zone in the phase space ( jWisdoml . 119871 ). Fig. 5 in (jWisdoml . 119871 ) shows 
that spatial rotation of Phobos with such initial conditions is not totally chaotic: there is a 
preferred orientation of the largest axis of satellite's figure in the direction to the planet. Let 
us consider an analogous graph for the chaotic motion of Prometheus close to s ynchron ous 
1:1 resonance. In the computation we adopt the values of b/a and c/b as given in EA'99l (see 
Table 1). The resulting projection of the spatial chaotic trajectory to the plane (0, 6 — f) is 
shown in Fig. |^a). 

Besides, in Fig. [9]^b), we build a three-dimensional density plot of the discrete projections 
of the trajectory to the plane (0, 9 — f). The output time step is taken equal to 0.01 of the 
orbital period. The square (0, 9 — f) G (— 7r/2,7r/2) x (— 7r/2,7r/2) is divided in a grid of 
40 x 40 pixels. The quantity N designates the number of the trajectory output points in a 
given pixel. 

The initial data is taken inside the chaotic layer of the resonance in the phase space 
section of planar rotation. The graph in Fig. |^a) is built for a trajectory on the time 
interval of 1 000 orbital periods, and that in Fig. Mjo) is built for the same trajectory on the 
time interval of 10 000 orbital periods. We present the graphs only for Prometheus; in the 
case of Pandora they look very similar. 

The semimajor axes of the orbits of Promethe us and Pandora are equal to 139 400 km 
and 141700 km, respectively f lFrench et all 120031 ): the mean radius of Saturn is equal to 
57 600 km. It follows then that the relative area of Saturn's disk as seen from the satellite 
(with respect to the area of the celestial hemisphere) is 8.9% for Prometheus and 8.6% for 
Pandora. These values give the average relative time that the largest axes of figures of these 
satellites would be oriented in the direction to Saturn, if orientations of the satellites during 
the chaotic "tumbling" were isotropic. For the chaotic trajectory presented in Fig. |9] we 
have calculated the values of the average relative time of orientation towards Saturn; they 
have turned out to be equal to ~ 30% for Prometheus and ~ 22% for Pandora for the time 
interval of integration of 10 000 orbital periods, i.e., the "isotropic norm" is exceeded 3.3 and 
2.6 times, respectively. From the plots in Figs. Efa) and (b) and these numerical estimates 
it is clear that there exists preferred orientation of the largest satellites' axes in the direction 
to Saturn, at least for the given test trajectory. 

More extensive additional test computations show that for the used time interval of 
integration of 10 000 orbital periods the values of the average relative time of orientation 
towards Saturn depend on the choice of initial data. The obtained values are in the range 
of 20-30%. The deviations may indicate that the observed anisotropy is a temporary effect 
due to specific initial conditions, and long term diffusion leads to its disappearance. This 
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b) 

Figure 9: Orientation of Prometheus in chaotic rotation, (a) Projection of the chaotic 
trajectory to the ((f), 6 — f) plane. The integration time is 1 000 orbital periods, (b) A 
three-dimensional density plot of the discrete projections of the trajectory to the plane (0, 
9 — f). The integration time is 10 000 orbital periods. 



remains an open problem. 



5 Conclusions 

We have studied possible rotation states of two small moons of Saturn, Prometheus and 
Pandora. There are two different regimes of synchronous rotation in the phase space of 
planar rotational motion of Prometheus: a-resonance and (3- resonance, i.e., it is subject to 
the "Amalthea effect" . Pandora has a-resonance only. Our analysis of stability of planar 
rotation of these satellites with respect to tilting the axis of rotation has shown that a- 
resonance for Prometheus is unstable or close to instability, i.e., the satellite most probably 
cannot reside in the given regime of synchronous rotation. Rotation of Prometheus in /3- 
resonance, as well as rotation of Pandora in its only possible a-resonance, is close to the 
attitude instability. Both satellites also possess period-doubling bifurcation mode of en- 
resonance in the phase space of rotation. Rotation of both Prometheus and Pandora in 
this mode is attitude unstable. With respect to multiplicity of synchronous states in the 
phase space, Prometheus is unique amongst the satellites with known inertial and orbital 
parameters. So, whether it is in chaotic rotation or not, its potential rotational dynamics 
are rich and complicated. To a less extent the same is true for Pandora. 

Our analysis of the attitude stability of planar rotation of Prometheus and Pandora for 
trajectories of various sort (periodic, quasiperiodic, chaotic) on a representative set of initial 
data, carried out by means of computation of the Lyapunov spectra, has shown presence of 
alternating concentric ring-like zones of stable and unstable trajectories around the major 
synchronous states for both satellites. Hypothetically, the belts of attitude instability might 
form "barriers" for capturing the satellites in synchronous rotation. 

In a numerical experiment we have demonstrated that the satellites in chaotic rotation 
can mimic ordinary regular synchronous behaviour: they can have preferred orientation for 
long periods of time, the largest axis of satellite's figure being directed approximately towards 
the planet. The presence of such anisotropy of orientation of the satellites in chaotic rotation 
might prevent clearing up the character of rotation in observations. Whether this anisotropy 
is a temporary effect due to specific initial conditions, with long term diffusion leading to its 
disappearance, remains an open problem. 
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